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Abstract 

The article is a translation of authors paper [I] printed earlier in the in- 
accessible edition and summarizing the results of research of gravitational 
waves damping problem in the cosmologic plasma due to the different 
interactions of elementary particles. 

Introduction 

The evolution of the cosmological gravitational waves (GW) in the Universe, in 
case that perfect fluid or a coUisionless gas is considered as a material medium 
model, is researched sufficiently [2] - 0- Let us note basic results: 

- perfect fluid doesn't change a vacuum character of the propagation of the 
gravitational radiation; 

- dispersion is essential only for super-density objects. 

There were made an attempts to provide an analysis of evolution for dissipative 
mediums (see, for example [I],[S]). In authors works [B], [5] — [TT] the damping 
of gravitational radiation during the propagation in the coUisionless gas was 
considered in the context of kinetic approach. Given paper may be considered 
as a concluding in the whole cycle of researches. 

1 Self-consistent model of weak gravitational waves 
evolution in the Friedmann Universe 

In the spatially flat Friedmann world with a metric: 

= a'{ri)[dTf - {dx^Y - {dx^f - [dx^f] (1) 
gravitational waves are described by metric perturbation: 

h^k = 9ik- 9ik= ^ifeWe-""'^"; a,/3 = T73, (2) 
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where = const - wave vector it's reference components, measuring by 

synchronic observer in metric arc ([T]), 



'a a 



where n — ^ n\ + n\ + ri^. It is obvious that in the flat Friedmann world waves 
of any length A can be counted like short: 

A < i? ^ cxD (3) 

where R — the radius of background curvature. Then on the hik additional 
conditions 2 of graduation and tracelessity can be imposed: 

/i,4 = 0, h = Q, h% = ^ h^np ^ 0; (4) 

o 

here and further the indexes lift and lower by the background metric 5ife, 
h = gikh^^ . How it will be seen from the further, the perturbation of the 
momentum-energy tensor (MET) of medium, determining by the GW appear- 
ance is representable in the form: 

STik = T{r])hik.. (5) 

Then the system of linearized by h Einstein equations can be written in the 
form 



where in further we will set: 

haisiv) = a{il)Sai3ip{ri), (6) 
at that Safi — const. Taking into account the background Einstein equations 
' a'\^ 8nea^ a" 47r 2/ „ \ ^tt 2 



, . , a'ie + 3p)~—a'e, (7) 

a J 6 a 6 3 

finally we will receive the evolution of GW conform amplitude equation ip{ri): 



-i(e-3p) + 2p + 2T 
6 



= 0, (8) 



where for convenience we have proceed to conserving on the ultrarelativistic 
stage conform densities of energy and pressure e = a'^{ri)e, p = a{ri)'^p etc. 
For calculation of t^t]) we will advert to relativistic kinetic equations 

[HJ] = J[f]. (9) 

The appearance of GW leads to the distortion of mass surface, thereof it is 
necessary to produce a renormalization of momentum [3]-[l]. It is simpler to 



^Here and further the Greek letters run the values from 1 to 3; Latin - from 1 to 4. 
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make an indicated renormalization by transformation Q of momentum vari- 
ables Pi to the momentum variables 

^^=P^- ^h'^pk, K = + ^/ifPfe + 0{h^). (10) 

The Jacobian of this transformation in case of traceless perturbations of metric 
in linear by h approach is equal to one. In new variables the Hamilton function 
coincide with the unperturbed value accurate to members 0{h^): 

H{x,p) = iff^W = \ 5^' P#fc =H {x,p). (11) 

Therefore g =9 +0{h ), differential of volume of momentum space is invariant 
with respect to transformation ([TT 



, 2s + l a^p 1 2^ 

2s + l a^p o 1 

=S{H --m)^dTT . (12) 



(27r)3 ^ 

Taking into account the fact, that momentum variables can be contained in 
invariant scattering amplitude Mif by means of all possible contractions of 
type {p,p'), it can be strictly shown, that collision integral accurate to 0{h^) is 
invariant with respect to transformations (|10|) . Thus, representing a distribution 
function f(x,p) in form 

f{x\p,)^ foiv,n) + Sf{x\Vk), (13) 

where /o(?7,P4) - isotropic solution of kinetic equations ([9|), and integrating the 
kinetic equations © , we will obtain linearized kinetic equations for distributions 
deviations Sf: 

^'^ + l^^^^^"^''-j'^fo,Sf]. (14) 

The isotropic part of distribution function /o satisfies to zero approximation of 
equation 

p°^ + -(PoP)|;^ = j"[/o]. (15) 

or] a o¥o 

Calculating the MET perturbation with the account of transformation pU)). 
we'U find: 

6T,k = -ph,k + P^^kSfdn, (16) 

where summation is carried out by all kinds of particles, participating in re- 
actions. From linear equations (IT4)) is clear, that Sf ^ S'c^P^P''. Actually, 
supposing such dependence we arrive at conclusion that collision integral in 
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consequence of it's invariance can be only a linear combination of contractions 
of type S'„;3F"P'^, 5'a/3n"P'^, Sap-nf^n^ , S^pS"'^ , from which only the first is dif- 
ferent from zero. And integrals of type (|16p are defined by the only selected 
spatial direction, n^. Therefore these integrals become the linear combination 
of contractions: 



from which only the last is different from zero and equal to Sap- 
Thus, finally we'll obtain: 



6Tik = -h,kp + hikTf, 



(17) 



where 



-y 



and the designation = SapSap is incorporated, - summation is carried out 
by repeating indexes. 

Thus, equation ^ takes form: 



--(e-3j5) + 2r/ 



= 0. 



(18) 



2 Calculation of collision integral 

Let's turn to calculation of r/. In case of equilibrium distributions /o an integral 
can be represented in form |llj : 

^M/o'5/] = ^4^K(P,r?), (19) 
i ± /o 

where 

K(P,77) = 

J] d7r5(P; - ¥,)W,f n /o 11(1 ± /o)./o)' (20) 

/,» / / 

Here summation is carrying out by all initial, (i), and final, (f), states of given sort 
particles, participating in reactions; ^ dir means the product of momentum 
volumes of all particles, except given, Wif - invariant scattering matrix (details 
see in Ref. [TB]). 

Let's illustrate a statement ([T5)) for four-particle reactions ah < — > cd. In 
that case collision integral takes form: 

1 



2(25a + l)(25'b + l) 

^^^(2^)45(p, + n - p, - f,)Jm;j\\ 
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X {fMil ± ± /,) - fahil ± /c)(l ± fd)} . (21) 

Multipliers (25a + (SS'i + 1)^^ (ED) correspond to the average by polarized 
states of particles [14 , and coefficients 1/2 respond to the selected normalization 
of colliding particles wave functions 15, . Substituting f — fo + 6f, where /o — 
isotropic equilibrium distribution function, and 6f ^ S'q/jP^P^, and linearizing 
an integral (I^TI) by 6f, it is easy to make sure that addends in the figure brackets 
of type Sfcf^{l ± /a)(l ± fb) during integration by dn turn to zero. Actually, 
directing an axis along the wave vector Ua we'll achieve, that Sfc will be 
proportional to expressions (P^2 ~ ^cs) ^c2^c3- Since multipliers at Sfc are 
invariant relatively to transformations P2 — > P3 and P3 — ^ P2 (for all parti- 
cles simultaneously), and Sfc at that change a sign, then during integration by 
momentums in infinite limits, this added turns to zero. Non-zero contribution 
in will give only members, containing Sfa, since integration doesn't carry 
out by momentums Pq. Thus, discarding non-sufficient for hot model statistical 
multipliers, we'll find for ultrarelativistic particles: 

K(ry,p) -i.(p,77)p, (22) 

where effective frequency of collisions ^'(p, 77) is equal: 

iy{p, 77) = 

00 4pij 

{2Sc + l){2Sd + l) 



l6TT^{2Sa + l)p2 



J fo{q)dq j satot{s)ds] (23) 




p - absolute value of physical momentum, 
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- total scattering cross-section, 

X^-t/s, S = {Pa+PbT, t={pa-PcY 

- kinematic invariants. 

When all interacting particles are ultrarelativistic, kinetic equations (|14l) 
take form: 

dSf Fc^dsf 

— h avoj — 



1 O/o So^pPcPp ff\' 

2 dP P V a 



dr] P dx" 

(24) 



It should be pointed out, that here P = Pq differs from physical momentum by 
multiplier (P = a{r])pphys)- It is obvious, that at > the second member 
in equation, conditioned by interparticle collisions, leads to the decrease of Sf, 
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i.e., to the relaxation towards the equilibrium distribution. Let's write down the 
solution of equation (l24l) . turning to zero at SajS — and owning the structure 
of coUisionless equation at 770 ^ 0: 



<5/ 



1 <SQ^]PJP^^_j^„a;°_i7r„ -7(,,) 



2 ap 


P 


X hm < 









(25) 



where unit vector tTq 



and damping decrement are incorporated: 



7(77, P)= iy{r],P)a{T])dr] ^ iy{t,P)dt 



(26) 



no 



It is not too hard to calculate the MET perturbation, conditioned by 6f, 



ttS, 



'^0 ^—iriaX^ 



4a2 



E 



2Sa + l 

(27r)3 



dP 



X lim { I — 

ijo-s-o V a 



+ J e-^^"'^^) (^^^ J[nir, - r^')]drj' 



(27) 



'70 



The summation in (j27p is carried out by sorts of particles a, participating in 
reactions, and incorporated function: 



sinx / 3 



J{x) = - 



having the asymptotics 
Ssinx 



3 cos a; 



J{x): 



, J(0) 



16 
15' 



(28) 
(29) 



Let's consider now the relationships, that are correct on the ultrarelativistic 
stage of Universe evolution (see., for example, [TB]): 



a = aif], 
^ {aT f 



45 2 



N - the statistical factor of particles number: 
iV = Es(25+l) + 1^^(25+1). 



(30) 
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Then we'll obtain: 



45 



OC 

y{2S+l) f dzz^^e'^C-^"); 
^-^ J oz 



X lim <; ( ^ ) J[n{j^ - r;o)] + 



(31) 



where z = V/T = p/T. 



Thus, finally we have - the evolution of cosmological GW in isotropic ultra- 
relativistic gas (e = ip) defines by equation: 

iP" + 'n^^P + 16na-^TfiP ^0, (32) 

where Tf is described by expression 1131]) . 

3 Extreme cases 

Let's at first consider the case of fast relaxation 

7 > 1. (33) 

Producing the asymptotical estimation of integral (1311) by Fourier method, we'll 
reduce an equation (|32]) to the form 



r+v'^+lr-) ^— = 0, (34) 

5 \rij aveffv 

where 

00 

Y^{2S + l)fdz z'^^, 
^ J dz v{z,7^) 

(35) 



4^(25+1) j dz z^fo 



In consequence of condition ([33]) the last member in the left part of equation 
P4p is small in comparison with the first two members - it responds to the weak 
damping of GW vacuum oscillations: 

" / T—2 ^ 
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i{nri- 



(36) 



At i^efft oo the damping of GW vanishes: in the ultrarelativistic fluid GW 
propagate as well as in vacuum. This fact was mentioned in the beginning of 
the article. 

Let now 7^1. In the coUisionless approximation an equation (|3ip can be 
simplified: 

mna^Tf = ^ lim i ( - ) J[r]{ri - 7/0)] + 



(37) 



This expression have two asymptotics. In the long-wavelength limit {nrj <C 1) 
we'h find 



Ibira Tf ~ - — 
5 



(38) 



Substituting (|38p in the equation |32l we'll obtain researched in [5] oscillations, 
as a solution: 



r PV3, ^ 
G+ cos In 7177 



^ . /3V3. ^ 
C_ sm in 717; 



(39) 



Such behavior, however, is forming only in the case that coUisionless situation 
occurred from the very beginning at 77 = 0. And if strong-collision phase pre- 
ceded the coUisionless phase (down to 77 ~ 7/0), then instead of ([55]) we have 
now: 
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V'(??) V'(?7o) 



T VoV 
and equation (j32[) becomes inhomogeneous 



8 V 



' V'('?o) 



5 7;2 5 7^977 



(40) 



(41) 



It's solution, satisfying to the sewing condition at the moment rj = rjQ^ has a 
form: 



^(^)=^(^o)^+^ /^sinj^ln^ 
Vo V '?o V 2 770 



(42) 



Main part of GW amplitude at that remains constant. Thus, a solution for 
long GW's, obtained in Ref. [5], is true only in the case that Universe started 
from the coUisionless phase, and is false in the case of existing of initial hydro- 
dynamic stage. At the last case GW stores it's final condition on hydrodynamic 
stage. 
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4 The evolution of short waves 

Let's consider the evolution of short {nrj ^ 1) GW's at arbitrary 7(77, P). In 
this case it is convenient to lay: 

iP = ^{Tjf exp (-i I ndrj) , (43) 



where fl{ri) - is a large value, ip{r]) - slowly changing function. At 77 — > GW's 
with any 77 are long, therefore it is necessary to redefine the solution (1251) at a 
point of time 770, when GW's become short: 770 > Carrying out necessary 
calculations, we'll find in the case of ultrarelativistic particles: 

00 

3 n ^2s + l /", 4^/, 



(44) 



loTra Tf — > ^ / azz < —n — 



f7' 



n 



where fl' = Q + ii>{r],p)a{ri). For the weak-collision plasma « 77, holding in 
Ref. pi)) members — 7^^, we'll reduce an equation ([5^ to the form: 

i^" - 2in^' - (r!^ - n^)^j - a'-^^I = 0, (45) 

7777-^ 

where v^ff defines analogically by (j35]). In mentioned approach the solution of 
equation ([55]) is: 

V'c.A;i/2(,^)^ exp (-170(7?)), ^ ^ 



where 



Let's point an inquisitive fact: neglecting damping of GW and introducing value 
$ = ip/a, we'll obtain from pS)) an equation for the mass scalar field 



- + 777q$ = 0, 



R 
"6 

with the effective mass of graviton 

1 

Let's clarify the question about GW evolution, which make sense only in the ap- 
proximation of geometrical optics. Averaged pseudotensor of momentum energy 
is equal to 

O/TT 
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Substituting here we'll obtain MET of perfect ultrarelativistic fluid with 
energy density 

2 2 

ec = exp[-7G(?7)], (49) 
i.e., the value (or k{r])) can be named as a decrement of GW energy damping. 

5 An example of GW energy damping decre- 
ment calculation in the reactions of type ee+ — > 
hadrons 

For illustration we'll carry out a calculation of GW damping according to stated 
above scheme for reactions of type ee^ — > hadrons, the cross-sections of which 
have a scaling behavior [17] 

-m(^) = ^E^N^^, (50) 

where e^- fundamental charges; their number, and, therefore the value Q, weakly 
depend from s. After calculations we'll find 

^eff{s) = TQ (T), (51) 

7a(oo) = ^(^^^ ^ 1. (52) 

As it was mentioned above, the GW damping is possible during synchronous 
fulfillment of conditions nrj > \ (short GW's) and ee//a < n (weak-collision 
gas). In the opposite case we have either long GW's, or GW's in the perfect 
liquid; in the both cases damping is absent. Therefore, regarding the fulfillment 
of these conditions in the initial point of time ryo, we arrive to the conclusion, 
that 7g(oo) ~ 4/3^, where /3 < 1 - unknown factor, guaranteing the severity of 
made approximations, moreover all GW's, having the wavelength 

on the actual point of time, damp with the same decrement (where - wave- 
length of relict photons). 

6 The results of GW damping decrement calcu- 
lation in the different reactions 

Let's diagram the evolution of GW in hot Universe as a dependence from the 
behavior of total cross-section of scattering. Let's isolate following typical situ- 
ations: 
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A: 


ni] ^ 1, 


avr/ ^ 1 


- long waves in the liquid, 


B: 


nr] 3> 1, 


avr/ ^ 1 


- short waves in the hquid, 


C: 


nij <^ 1, 


avr/ < 1 


- long waves in the weak-collision gas, 


D: 


nij » 1, 


avrj < 1 


- short waves in the weak-collision gas 



riQ 



n> no C ^ (b) 



n < no C ^ A 
<m C ^ (5) 



B 



^tot 



^tat 

no ~ ■ 



> m A — > B 



n> no A^ B @ 



n < no A 



C 



Situations, where coUisional GW damping occur are circled; the domain of GW 
amplitude fixation are squared. The problem of damping effectiveness requires 
the realization of calculations for concrete reactions. 

The results of such calculations, fulfilled in papers [9]- [11], are represented 
in the Table 1. In two last columns of the Table represented the results, sum- 
marized by all reactions of given type. G - Fermi constant, 

0, laxN'^nipi 

o- = 

nix 

nipi - Planck mass, mx - mass oi X - bosons, ax - interaction constant, Nx - 
the number of particles participating in reactions. 

Summing the results, we'll note, that collision damping of cosmological GW's 
in the all conceivable reactions is not great. However, the reason of weakness 
of GW damping isn't the weakness of interparticle interaction, but the circum- 
stance, that exactly in the medium with intense interparticle interactions, i.e. 
in the perfect fiuid, GW's don't damp. In order to provide an appreciable GW 
damping, collisions don't need to be very frequent, that, in turn, leads to the 
weakness of damping. 
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Table 1. The damping decrement of cosmological gravitational waves 
in the early hot universe for different processes in cosmological plasma 



Process 



Total cross-section of 
scattering 

crtotis) 



The region of mod- 
ern wavelengths, where 
damping is maximal 
(in cm) 



Maximal value 70 



ee"^ ^ 77 



In- 



100 



0,01 



ee^ ^ fifi^ 



~37 



4/32 



ev O ev 



3) 



ev o eiy 



G%s/3tt 



2 • 10 



20 



0,74 



|M|2 = 87ra:,m2iV 

CT > 1 



gV3 



0,29^ - 3 • 10^2 



i)By data of Ref.[T2]; 2)by data of Ref . [9] , [10] ; 3)by data of Ref. [12]; 4)by data 

ofRef. \m,\m- 



The unique sufhciently effective (in order of GW damping) interactions are 
the interactions with scaling behavior of cross-section. Exactly scale-invariant 
interactions in consequence of identical time law of collision frequency with GW 
frequency can infinitely long influence on GW and thereby apply to the sufficient 
damping of GW. However at attempt to calculate a GW damping decrement in 
the most effective region v ^ k we'll fall into long-length collisionless phase, in 
which the GW energy determination is difficult to define unambiguously. The 
calculations of damping decrement during electro-weak interactions generally 
confirmed the Hawking's estimation [7] and rather improve it. Realized calcu- 
lations point that the height of energetic spectrum of cosmological GW's can 
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be decreased approximately in 1,5 - 2 times at lengths A > W^P^ /N'Q^ cm and 
additional to this the same effect at A ~ 2 x 10^°cm. 
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